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SOME  THEORETICAL  ASPECTS  OF  FIRE  INDUCED  FLOWS  THROUGH  DOORWAYS 
IN  A ROOM-CORRIDOR  SCALE  MODEL 


James  G.  Quintiere  and  Karen  DenBraven 
Abstract 

Fire  induced  flows  were  measured  in  a 1/7 th 
scale  model  room-corridor  by  measuring  velocity 
and  temperature  profiles  in  the  room  and  corridor 
exit  doorways.  The  corridor  exit  door  width  was 
varied  by  a factor  of  10  as  an  experimental  pa- 
rameter while  the  average  fire  room  temperature 
was  held  constant.  The  applicability  of  an  ori- 
fice flow  model  to  predict  the  doorway  flows  was 
examined.  The  ratio  of  measured  to  theoretical 
mass  flow  rate,  defined  as  the  flow  coefficient, 
was  found  to  vary  with  door  width  and  flow  direc- 
tion. The  coefficient  ranged  from  0.4  to  slightly 
greater  than  1.  No  explanation  has  been  developed 
to  account  for  this  variation. 

Key  words:  Corridor;  experiment;  room;  scale-model; 

theory;  compartment  fires;  buoyancy  flows;  flow 
coefficients. 


1.  INTRODUCTION 

It  has  been  found  fruitful  to  study  fire  induced  flows  in  a scale 
model  facility.  The  model  consisted  of  two  compartments  separated  from 
each  other  and  the  surroundings  by  two  doorways  - representing  a room 
and  an  adjoining  corridor  (figure  1).  The  first  study  with  this  facility 
indicated  the  occurrence  of  four  alternating  horlzontol  flow  layers  in 
the  corridor  [1].  A subsequent  study  by  McCaffrey  and  Quintiere  [2] 
illustrated  complex  three  dimensional  flows  in  the  corridor,  and 
demonstrated  that  the  multi-layer  corridor  flow  pattern  also  occurs  in 
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a full-scale  counterpart  similar  to  the  model.  A more  recent  study  [3] 
further  illustrated  these  corridor  flow  patterns  through  visualization 
techniques,  and  demonstrated  some  limitations  of  theoretical  flow  models 
in  coping  with  the  prediction  of  these  flows. 

Several  theoretical  approaches  have  been  made  to  predict  fire 
induced  flows  between  compartments.  One  approach  relies  on  solution  of 
the  field  equations  and  a suitable  turbulent  model  [4].  Another  approach 
relies  on  a model  that  depicts  the  flow  into  an  upper  hot  layer  of 
uniform  temperature  with  a lower  cold  layer  [5].  The  doorway  flows  are 
treated  as  orifice-like  flows  as  described  by  Prahl  and  Emmons  [6]. 

The  latter  approach  of  reference  [5]  requires  a knowledge  of  doorway 
flow  coefficients,  and  a quantitative  prediction  of  flow  rates  entrained 
into  the  jets  of  flow  issuing  through  the  doorways.  Neither  approach 
has  been  sufficiently  developed  to  permit  appropriate  comparison  with 
experimental  measurements. 

This  paper  examines  the  applicability  of  the  orifice-flow  two-layer 
fluid  concept  for  doorway  flows  measured  in  the  room-corridor  scale 
model.  The  independent  variable  in  this  study  was  the  width  of  the 
corridor  exit  doorway.  Consequently,  the  flow  characteristics  will  be 
expressed  as  a function  of  the  exit  doorway  width.  The  theoretical 
model  for  the  doorway  flow  rates  is  developed  as  a function  of  the 
temperature  distributions,  from  which  the  flow  coefficients  are  determined. 

2.  EXPERIMENTAL  DESIGN  AND  ANALYSIS 

Experiments  were  done  in  a small  (1/7)  scale  model  corridor.  The 
model  was  0.35  m wide  by  0.35  m high  by  1.3  m long,  with  a variable 
width  (corridor  exit)  door  opening  at  one  end  (figure  2).  Within  the 
model,  a 0.35  m long  burn  room  was  set  up  at  the  end  opposite  the  exit 
doorway.  A 0.076  m by  0.140  m porous  plate  diffusion  flame  gas  burner 
was  placed  on  the  centerline  along  the  wall  opposite  the  room  entrance 
to  simulate  a room  fire.  A small  pitot-static  tube  was  placed  in  each 
of  the  doorways  for  velocity  measurement  along  a vertical  traverse. 


2 


0.14  m from  the  side  wall,  with  thermocouples  attached  to  record 
temperatures.  All  thermocouples  were  made  from  fine  5 mil  (0.127  mm 
diameter)  chromel-alumel  wire,  and  were  aligned  horizontally  facing  into 
the  flow.  Thermocouples  were  also  placed  within  the  room  and  the  corridor, 
0.04  m from  the  sides. 

The  corridor  model  was  supported  above  the  laboratory  floor  so  that 
the  exit  doorway  did  not  have  a floor  that  extended  beyond  the  corridor. 
However,  a floor  was  added  for  some  cases  with  no  significant  change  in 
the  results.  The  most  significant  effect  of  the  surroundings  on  the 
results  was  caused  by  the  laboratory  ventilation  conditions.  At  times 
"windy"  room  conditions  prevailed  which  necessitated  terminating  the 
run.  On  those  occasions,  the  laboratory  ventilation  system  was  shut 
down  resulting  in  relatively  "still"  room  conditions,  and  a new  run  was 
initiated . 

During  an  experiment,  the  burner  was  lit,  and  the  room  and  corridor 
were  allowed  to  heat  until  the  ceiling  layer  in  the  room  reached  a 
stable  temperature  of  330  +30  °C.  (Maintaining  a constant  room  tempera- 
ture tends  to  minimize  the  effect  of  temperature  on  flow.)  Vertical 
traverses  were  then  made  in  the  doorways,  room,  and  corridor  with  the 
pitot  tubes  and  thermocouples. 

Data  were  recorded  for  the  inflow  in  Increments  of  0.013  m from  the 
corridor  floor  up  to  the  neutral  plane  as  indicated  by  a zero  reading  by 
the  pitot  tube,  i.e.  zero  velocity.  These  measurements  were  only  taken 
on  or  near  the  doorway  centerline.  Data  for  the  outflow  were  recorded 
in  a similar  manner  using  a pitot  tube  and  thermocouples.  Centerline 
measurements  were  taken  in  increments  of  0.013  m,  and  in  addition 
measurements  off  the  centerline  were  taken  every  0.025  m (figures  3 and 
6).  Exit  doorway  opening  width  was  varied  from  0.034-0.35  m for  both 
inflow  and  outflow,  while  the  room  doorway  height  and  width  remained 
constant . 
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To  obtain  estimates  of  the  orifice  coefficient  C,  the  ratio  of 
experimental  to  theoretical  mass  flow  rates  were  calculated.  Experimental 
velocities  and  mass  flow  rates  were  obtained  from  the  doorway  temperature 
and  pitot  tube  data.  Theoretical  velocities  and  mass  flow  rates  were 
calculated  from  measured  temperatures  as  indicated  in  the  following 
analysis . 


2.1.  Doorway  Flow  into  a Compartment 

Bernoulli's  equation  was  applied  to  a streamline  starting  from  rest 
outside  the  dooi-way  to  a position  downstream  of  the  doorway  where  the 
streamlines  are  parallel  (this  can  be  considered  a vena  contracta) . 

This  is  illustrated  in  figures  4 and  5,  where  the  velocity  at  the  vena 
contracta  of  the  jet  is  given  by 


Vj (y)  = 


P^(y) 


(1) 


here  is  the  density  of  the  incoming  fluid.  The  pressure  in  the  jet 
is  equal  to  the  pressure  p^  in  the  stagnant  fluid  surrounding  the  jet. 
C^,  the  viscous  coefficient,  is  close  to  1 for  pipe  flows. 


Then  taking  the  pressure  distribution  to  be  hydrostatic,  p^,  the 
corresponding  pressure  outside  the  doorway  at  y for  fluid  with  no 
motion  or,  alternately,  the  stagnation  pressure  in  the  doorway  at 
position  y,  can  be  expressed  as 


y 

Pq(y)  = + g J pQ(y)dy  (2) 

o 

and  p^,  the  static  pressure  just  inside  the  doorway  at  y associated 
with  the  static  fluid  surrounding  the  vena  contracta  of  the  incoming 
jet,  is  given  as 


y 

Pi(y)  = Pj^  + g f P^(y)dy 
o 


(3) 
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where  y increases  from  the  neutral  plane  down,  and  is  the  pressure 
at  the  neutral  plane,  (the  height  at  which  p^  = p^) . 

is  the  density  of  the  stagnant  fluid  surrounding  the  jet. 

Analogous  to  the  development  of  orifice  flow  theory,  the  jet  cross 

sectional  area  A,  at  the  vena  contracts  is  related  to  the  doorway  inflow 

area  A by  a contraction  coefficient  C , i.e. 

c 

A.  = C A (4) 

J c 

The  assumption  was  made  that  the  contraction  coefficient  is  constant  and 
uniform  over  the  flow  field  for  each  streamtube  between  the  door  and 
the  jet.  By  conservation  of  mass  for  each  horizontal  layer  of  stream- 
tubes,  this  leads  to 

v(y)  = v^(y)  (5) 


The  accuracy  of  this  assumption  will  be  examined  by  comparing  the  theo- 
retical results  for  v with  measured  values. 


Combining  equations  (1)  and  (5)  yields 


v(y) 


(Pq(y)  - P^(y))dy 

PoW 


(6) 


where  C = C C . 

V c 

Assuming  a perfect  gas  for  air. 


PT 


= P 


T 

00  00 


(7) 


where  p and  T are  reference  state  values. 

00  oo 
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Substitution  in  equation  (6)  yields 


or 


(8) 


v(y)  cJ  2g  T^(y)  J T.  (y)^ 


(9) 


Since  the  measurements  were  in  the  form  of  T as  a function  of  the 
distance  z from  the  floor  a change  of  variables  was  made  such  that 


dz 


(10) 


where  z is  measured  from  the  floor  and  z is  the  height  of  the  neutral 

n 

plane. 


For  this  experiment,  was  the  temperature  in  the  doorway  entrance 
to  the  room  or  corridor,  and  was  the  temperature  in  the  corner  of  the 
appropriate  compartment.  Therefore,  to  calculate  flows  into  the  corridor 
exit,  and  were  used.  Simpson's  1/3  rule  was  applied  to  carry  out 
the  required  numerical  integration.  This  method  requires  an  odd  number 
of  data  points.  For  locations  along  the  vertical  traverse  where  only  an 
even  number  of  data  points  were  available,  the  last  two  values  at  the 
upper  limit  of  the  integral  were  averaged  to  create  another  data  point 
between  the  two  which  was  then  used  to  complete  the  calculation.  The 
mass  flow  rate  was  determined  by  integrating  the  mass  centerline  velocity 
pv  over  the  height  of  the  inflow.  The  experimental  mass  flow  is 


m 

exp 


p (z)  V (z)  dz 
o exp 


(11) 


o 


where  W is  the  doorway  width.  Likewise  the  theoretical  mass  flow  is 
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From  these  mass 
was  calculated. 


o 

flow  calculations, 


p (z)  v(z)  dz 
° “C~ 

the  flow  coefficient 


for 


(12) 


the  inflow 


C 


m 


in 


exp 

th 


2.2.  Doorway  Flow  from  a Compartment 


(13) 


The  theoretical  velocities  in  the  outflow  were  calculated  by  assuming 
that  the  flow  just  outside  of  the  opening  is  of  a uniform  pressure  over 
a horizontal  slice  and  equal  to  the  surrounding  pressure  p^  at  the  same 
height  (figure  6).  From  Bernoulli's  equation 


where 


Pi 


is  the  stagnation  or  upstream  pressure  at  height  y inside 
the  compartment 


( )j  denotes  the  position  in  the  jet  emerging  from  the  opening 
where  p . = p . 

J o 


Again,  if  the  neutral  plane  has  pressure  p^,  and  assuming  horizontal 
velocity  only  so  that  the  vertical  pressure  variation  is  hydrostatic,  then 


dy 


(15) 


o 

where  y is  the  distance  increasing  upward  from  the  neutral  plane,  and 


so  that 


(16) 
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(17) 


Pi  - Po  = 8 


/ ^1=0  " 


Since  p.  = (the  density  in  the  compartment  just  upstream  of  the 
doorway)  the  combination  of  equations  (14)  and  (17)  gives 


(18) 


Also,  if  there  is  a contraction  in  the  jet,  as  in  orifice  pipe  flow,  then 

A.  = C A.  (19) 

J c 

If  this  contraction  coefficient  is  uniform  over  the  exit  jet  then 


V = C V . , 
c 3 

and  if  C = , then  equations  (18)  and  (20)  give 


(20) 


(21) 


Or  since  p . = p T /T . 

1 o o 1 


v(y)  = C 


^ 2g  T,<y)  jf  ^ 


which  can  be  compared  to  equation  (9)  for  inflow. 


dy 


(22) 


As  in  the  case  of  inflow,  experimental  velocities  were  obtained 
from  pitot  tube  measurements.  To  calculate  the  doorway  flow  coefficient 
C,  experimental  and  theoretical  mass  flows  were  again  obtained. 

H-z 

1 

v(y)  dy  (23) 


*th  ■ ” 


/ - 


or 
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and 

W H-z 

I f P^(y)  v(y)  dy  dx  (25) 

o o 

Here,  the  integration  was  carried  out  over  the  entire  region  of  outflow, 
since  data  were  available  and  a uniform  mass  velocity  profile  across  the 
doorway  width  does  not  appear  to  be  an  acceptable  assumption.  An  example 
of  experimental  results  for  the  exit  outflow,  p ^v  , is  shown  in  figure  7. 
These  results  correspond  to  the  data  given  earlier  in  figure  *2.  Hence, 
from  the  above,  the  flow  coefficient  is  defined  as 

m 

C = ^ (26) 

"'th 

3.  RESULTS 

Figure  8 shows  the  results  for  doorway  mass  flow  rates  taken  in 
these  experiments  and  in  previous  experiments  with  the  corridor  model 
[1,2].  The  results  indicate  that  as  the  exit  width  is  increased  the 
exit  flow  increases  almost  directly,  yet  the  room  flow  increases  only 
slightly.  The  difference  between  the  doorway  flows  represents  the  net 
entrainment  or  recirculation  that  occurs  in  the  corridor.  No  adequate 
theory  has  been  developed  tp  predict  this  mixing  rate. 


Figure  9 shows  the  results  of  calculations  for  the  flow  coefficient 
C vs  Reynolds  number.  This  is  compared  with  the  experimental  results 
from  Prahl  and  Emmons  for  kerosene/water  doorway  flows  [6].  Results 
appear  to  agree  qualitatively  but  not  quantitatively.  The  Reynolds 
number  is  defined  as 

m 

Re  = (27) 
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where  y is  the  viscosity  at  the  average  temperature  of  the  fluid.  In 
this  form  the  characteristic  length  dimension  is  the  flow  height. 

Figure  10  shows  the  exit  doorway  coefficients  as  a function  of  the 
exit  doorway  width  ratio  (W^/W^) . These  results  show  the  reproducibility 
of  the  data  for  different  experiments,  and  the  variation  of  the  flow  co- 
efficient with  the  exit  doorway  width.  No  explanation  has  been  developed 
to  account  for  this  variation,  nor  to  explain  the  higher  flow  coefficient 
for  the  outflow  than  for  the  inflow.  Also,  although  not  shown  in  figure 
10,  the  room  doorway  at  a fixed  width  had  a flow  coefficient  for  the  flow 

into  the  room  of  about  0.77  for  all  values  of  W /W  . 

h u 

In  figures  11-14,  measured  centerline  temperature  and  velocity 
profiles,  and  the  results  of  the  calculation  of  theoretical  velocities 
are  presented  for  four  typical  experiments,  where  the  corridor  exit  door 
opening  width  (W  ) ranged  from  0.071  Am  to  0.35  Am.  The  data  were 

Cj 

taken  up  to  the  height  of  the  neutral  plane  z^.  From  the  temperature 
data,  it  can  be  seen  that,  as  the  exit  doorway  is  narrowed,  more  hot 
fluid  enters  the  room  doorway.  This  is  equivalent  to  imagining  that  the 
height  of  the  hot  upper  layer  of  fluid  in  the  corridor  falls  below  the 
neutral  plane  height  in  the  room  doorway.  The  neutral  plan-e  is  lower 
for  the  room  doorway  than  it  is  for  the  exit  doorway  in  all  cases. 

These  features  are  illustrated  in  figures  15  and  16.  The  height  of 
an  effective  hot  layer  was  determined  by  a fixed  temperature  (50  °C  for 
the  corridor  near  the  exit  doorway,  and  60  °C  for  the  flow  entering  the 
room)  or  by  a tangent-intercept  method  illustrated  in  figure  16  for  the 
exit  doorway.  In  general,  both  methods  appear  to  give  the  same  results. 
However,  the  results  for  the  thermal  layer  given  in  figure  15b  are 
probably  more  Indicative  of  the  entire  corridor.  In  fact,  these  results 
are  in  agreement  with  corridor  smoke  layer  heights  determined  in  an 
earlier  study  [3].  The  thermal  heights  indicated  for  the  room  doorway, 
except  for  small  W /W  cases,  are  probably  not  indicative  of  the  corridor. 
A thermal  layer  in  the  corridor  lower  than  the  room  doorway  neutral 
plane  could  result  for  large  W /W  from  conduction  between  the  slowly 

£j  D ^ ^ 


moving  hot  and  cold  flows.  However,  it  is  more  likely  to  be  an  in- 
correct result  due  to  radiant  heating  of  the  thermocouple  from  the  hot 
room  and  fire. 


Since  fire  induced  doorway  flow  models  incorporated  in  global  room 
fire  models  have  been  based  on  a conception  that  the  fluid  is  stratified 
into  two  zones  of  uniform  temperature,  it  was  of  interest  to  examine  our 
data  from  this  point  of  view.  The  corridor  doorway  inflow  was  determined 
using  this  two-layer  temperature  model.  The  equation  follows  from  equa- 
tions (10)  and  (12)  provided  T^  is  taken  as  being  equal  to  T^  up  to  an 
effective  corridor  thermal  height  z^.  The  result  is 


zero  velocity  position.  The  value  for  C was  equal  to  the  corresponding 
value  determined  from  the  velocity  and  temperature  data,  i.e.,  equation 
(13).  The  results  of  this  calculation  were  compared  to  the  measured 
values  of  flow  rate  (figure  17).  In  general,  the  results  are  in  good 
agreement,  but  the  value  calculated  from  equation  (28)  tends  to  under- 
estimate the  experimental  flow  rate  in  a few  cases  by  up  to  30%. 


It  is  difficult  to  generalize  from  the  results  of  this  study  since 
the  effect  of  scale  may  limit  its  application.  However,  the  following 
specific  conclusions  can  be  made. 

1.  The  doorway  flow  coefficient  appears  to  vary  systematically  with 
doorway  width,  but  no  explanation  has  been  developed  for  this 
effect.  Also,  the  coefficient  appears  to  be  consistently  greater 
for  the  outgoing  flow  than  the  inflow. 


The  value  for  T^  was  determined  as  indicated  in  figure  16  as  the  maximum 
temperature  of  the  fluid  inside  the  doorway.  The  height  z^  was  deter- 
mined by  the  tangent-intercept  method,  and  z^  was  determined  from  the 
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2.  The  theoretical  results  for  velocity  are  in  good  agreement  with 
the  experimental  measurements.  This  may  suggest  that  a constant 

flow  coefficient  can  be  applied  to  calculate  the  velocity  distribution. 

3.  Although  a thermal  layer  does  not  exist  as  a discontinuous  inter- 
face as  idealized  in  a two-layer  model,  it  can  be  defined  from  the 
data.  Calculations  of  mass  flow  rate  based  on  experimentally  deter- 
mined layer  heights  did  yield  results  within  30  percent  of  the 
measured  flows  for  the  exit  doorway  inflow  in  this  experiment. 
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Figure  ].  Pictorial  of  scale  model  facility 
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Figure  3a.  Temperature  data  for  the  flow  out  of  the  corridor 
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Figure  4.  Calculation  of  Inflow  velocities 


18 


Figure  5.  Diagram  of  vena  contracta  for  inflow  (horizontal  section) 
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Figure  6.  Calculation  of  outflow  velocities 
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Figure  7.  Mass  flow  rates  for  flow  out  the  exit  doorway 
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Figure  8.  Experimental  mass  flow  vs  corridor  exit  door  width  (W  /W, 
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Figure  9.  Orifice  flow  coefficient  vs  Reynolds  number 
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Figure  iO.  Orifice  flow  coefficient  for  exit  doorway  vs  corridor 
exit  door  width  (W  /W  ) for  inflow  and  outflow 
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Figure  12.  Experimental  temperature  and  velocity  results 
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Figure  13.  Experimental  temperature  and  velocity  results 

and  calculated  velocities  for  W„/W^  = 1.0 
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Figure  14.  Experimental  temperature  and  velocity  results 

and  calculated  velocities  for  W„/W^  =0.6 
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Neutral  plane  height  and  effective  height  of  hot 
layer  based  on  room  doorway  temperature  vs 
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Figure  15b.  Neutral  plane  height  and  effective  height  of  hot 
layer  based  on  corridor  temperature  adjacent  to 
exit  doorway  vs  W /W 
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Figure  16.  Example  of  determination  of  z^,  z , T,  and  T from 
temperature  data  ^ ° 
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Figure  17.  Calculated  mass  flow  rate  vs  experimental  mass 
flow  rate 
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